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Abstract 

•/^ I We consider a spin-^ chain with competing nearest and next-nearest neigh- 

o 

^ , 

, an equiavelent to the ANNNI model. When studing thermodynamics of the 



bor interactions within a transverse magnetic field, which is known to be 



J2 I 2D ANNNI model Villain and Bak arrived to a free fermion approximation 

, that neglects heavy excitations from the ferromagnetic ground state, which 

Q , is an appropriate description close to the paramagnetic-ferromagnetic transi- 

O 

tion. In the vicinity of the floating-phase/anti-phase transition another sort 

•i-H . 

^ I of quasiparticles, but free fermions too, appears to be convenient. Although 

free fermions are a suitable tool for investigation of the phase diagram and the 
critical properties, they are defined on the fictitious lattice which makes the 
analysis non-rigorous. Here we deal with a proper fermion scheme which is 
especially effective for performing exact diagonalization calculations for cyclic 
chains. Systems up to size L = 32 has been analysed and the predictions of 
the effective fermion Hamiltonian has been confirmed. Various predictions for 
the infinite system and the critical properties are derived. 
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I. INTRODUCTION 

Frustrated Ising models within a transverse field have been investigated for a long time0. 
The transverse field plays the role of a tunable parameter by which one can induce a so 
called "quantum phase transition" at zero temperature that is driven by quantum fluctuation 
alonei (as opposed to a conventional, thermally driven phase transition). Frustration can be 
introduced via disorder, as for example in the case of quantum Ising spin glasses that have 
gained much interest quite recently!. However, it can also be produced in a regular fashion 
as e.g. in the anisotropic next nearest neighbor Ising model (ANNNI)EI, where competing 
nearest and next-nearest neighbor interactions are the origin of the richness of the phase 
diagram. This is the model that we intend to re-investigate in this paper. To be concrete, 
we consider the Hamiltonian, 7i, of a Id spin-| chain in a transverse magnetic field, which 
consists of the two parts: the classical Id ANNNI model plus quantum fluctuations imposed 
by transverse field: 

Ti- = Hcl + i/qu, (1) 

i^ci = -JE^r<+i + '^^E^r<+2 and H^^ = -Tj2at (2) 

i i i 

r = corresponds to a classical ANNNI model in which J is supposed to be positive, cr^ 
and are the Pauli matrices: 

/, n \ /n ,\ 

(3) 

The most interesting region of this model is the region around k = 1/2. The classical part 
Hcl has an infinite ground state degeneracy at this point, separating the ferromagnetic (FM) 
region (k < 1/2), where the ground-state is given by "all-spins- up" or "all-spins-down", from 
the anti-phase region (/t > 1/2), where the ground-state has a period 4 with two down-spins 
following two up-spins. Afteri < 2 > is a traditional notation for this anti-phase. Switching 
on the quantum fluctuations via a nonvanishing transverse field induces the presence of 
different phases: fixing F ^ J and increasing k in the vicinity of 1/2, one can distinguish 







1 






1^ 


a' = 




















.1 





. FM («) ^ 0); 

• PM, i.e., paramagnetic phase (exponentially decaying spatial spin-spin correlations, 
no long-range order); 

• FP, i.e., floating phase (algebraically decaying spatial spin-spin correlations, accom- 
panied with a modulation continuously changed with n, no long-range order); 

. <2>((a,^<+2)<0, («+4)>0). 

The FM-PM transition is well understood: It is within the same universality class as the 
two-dimensional Ising model or Id unfrustrated Ising model in a transverse field (see, e.g.,i'i). 

Villain and Baki in their seminal work on the two-dimensional ANNNI model argued 
that the < 2 >-FP transition is expected to be of the Prokovsky-Talapov type0 and that 
the FP-PM transition is expected to be of the Kosterlitz-Thouless typeini. Since this work 
is based on various plausible but not rigorously proven assumptions many attempts have 
been done to check these predictions with Monte Carlo simulations or exact diagonalization 
studies (for a review seei). However, for principle reason that we also try to clarify in this 
paper, such an endevour turns out to be very difficult and fails to provide the theory either 
with a conclusive support or with a clear falsification. 

A principal question of this work is how to catch the essential physics in the vicinity of 
the FP-< 2 > and FP-PM transitions. For doing this we present a regular expansion of 
Hamiltonian (0) in powers T/J', which allows to perform exact numerical diagonalization 
for longer chains than it is possible to achieve by a straightforward procedure. 

The organization of the paper is as follows: In Section II we recapitulate the free fermion 
picture of Villain and Baki before we present the above mentioned effective Hamiltonian. 
Then in Section III we compare the analytical predictions of this theory for finite systems 
with the results of exact diagonalization studies and proceed to extract the desired informa- 
tion about the critical behavior. Section IV summarizes our results. 
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II. FREE FERMION PICTURE 



Near k—1/2 the elementary excitations can be distinguished as hght and heavy. The 
main idea is to map the initial Hamiltonian onto the states where the heavy excitations are 
excluded from. 

We start with a conditional FM vacuum state, for example, ..+ + + + +.., then the 

excitation with an isolated (— ) spin, e.g., .. + H h + .. costs the energy AJ{1 — k) which 

is not small as k — > 1/2. However, the excitation with two or more sequential spins rotated, 

e.g., .. + H h +.. costs the energy AJ{1 — 2k) which vanishes with k — 1/2 . At the 

next step we can introduce quasiparticles, which are the domain walls (DW's) defined on 
the dual lattice. The latter coincides in Id with the middles of the links. The energy of a 

dw 

single DW, i.e., .. + + + | .. , is determined by the classical part of Ti.: 

e = 2J{1 - 2k) . (4) 

The DW's which occupy the nearest sites of the dual lattice repel each other with the cost 
of energy V = AJ'k. In general any state now is characterized by positions of the DW's on 
the dual lattice 

{£1,4,..., 4}. (5) 

For convenience we set the coordinates on a dual lattice to be integer numbers, 1, . . . , L, 
while the sites of a real lattice run half- integer numbers, say, 1/2, ...,L — 1/2. The quantum 
part of H plays a role of the kinetic energy of quasiparticles. In fact, applying H^^ to the 
state with a DW located at the site i, the spin from either its right or its left is changed 
by sign, that means a shift of a DW by one unit. Applying Hq^ to the site with no DW's 
surrounding it creates a couple of DW's. Hermitian conjugation corresponds to annihilation 
of those DW's. 

All the matrix elements described in terms of DW variables on the dual lattice can be 
summarized in the following Hamiltonian: 
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n = Ho + Hu Ho = VY. 

j 

3 3 

where r's are usual Pauli matrices and = (1 - )/2. A conditional vacuum is a state 
with no DW's (rj = 1 or n^-^'' =0 in (|^)), creation (annihilation) of a DW is realized by r~ 
(r+). A standard derivation of Hamiltonian (|^) and transformation from spin-operators a's 
to spin-operators r's is given in Appendix A. 

A routine procedure of the Jordan- Wigner transformation 

Tj = Cj exp Uvr 2^ nfc 
\ fc=i 

allows to deal with fermionic variables. Below we shall use the periodic boundary conditions, 
that means that the L-l-l-th site should be identified with the first. Hamiltonian can be 
rewritten now as 

L-1 

1-1= [euj + Vrijiij+i - T{c]c]j^i + cj+iCj + c]cj+i + c]+iCj)) + 

+ VniTiL - T{c\c\ + ClCi - c\cl - c\ci) exp{ink) (7) 

where k = Y^i n-j^ a total number of DW's. Evidently, k should be even number on a cyclic 
chain that results in exp(mA;) = 1. 

Assuming (cj T) one can obtain the effective Hamiltonian, 7i^^ , which reflects low 
energy properties of 7i with the energy scale of order (e, F). Note, that the l^-term makes 
two DW's energetically unfavorable if they occupy the nearest sites. Simultaneously, terms 
r(r+r+ + T^T^) of (|^) (or terms T{c'c^ + cc) of ([^)) should be excluded from the effective 
Hamiltonian, which now reads in r-variables 

= E^^f - r(r;rr^i + rrr+ ,)) - r(r+r, + r^rt) (8) 
i=i 

and in fermionic variables 

L-1 

^eff = E - r(c]c,+i + c]+iC,)) + T{C\CL + cici). (9) 



The constraint should be imposed on a possible wave function: DW's, or fermions, occupying 
the nearest sites, are forbidden. 

Let us try a wave function of the form 



t 



r' r' 



(10) 



where k should be even. We may search for the amplitudes /'s in a form of Bethe substitu- 
tion: 



f {mi, 7712, ■ ■ ■,mk) = '^^pexpi{qp^mi + . . . + gp^mfc) 

{P} 



(11) 



where {P} is a permutation of numbers {1, 2, ... , k}. Using the results of Appendix B we 
obtain a general expression for the eigenfunctions of Hamiltonian (0) in the V ^ oo limit 



mfe 



gig2mi g«<?2(»n.2-l) 



^jgi (mfe-fc+1) 



3«(j2(mfe-fc+l) 



3«'?fe("T-fc-'i-'+l) 



gt gt gt 
mi m2 ■ ■ ■ mj; 



|0). 



(12) 



The wave function (|T2D can be interpreted as a wave function of a tight-binding fermion 
model on a fictitious lattice. The coordinate of a fermion in a fictitious lattice coincides 
with that one in a real lattice minus the number of fermions situated from its left. We shall 
also use another interpretation which will be convenient in numerical diagonalization. It is 
consistent with introducing two kind of "particles", A and B. A is composed of a DW (or 
fermion) with a nearest empty site from its right attached. B represents an empty site which 
has no a nearest DW from its left. The A — B representation will be discussed in Section 
III in detail. 

The ground state energy of any intermediate state on the phase diagram between the 
FM and <2> boundaries 



E, 



ke-AVY, 



^ (2m-l)7r 



cos ■ 



m=l 



L-k 



ke-2T 



simr k/{L — k) 
sinyr/ {L — k) 



(13) 
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is a function of k, the total amount of DW's, which characterizes modulation of a spin 
structure. Eq.(|l^) describes the ground-state energies of the FM structure (fc = 0,Egs = 0) 
and of the <2> structure {k = L/2,Egs = Le/2) as well. Both structures, FM and <2>, 
can be unstable with respect to formation either DW's {k = 2) or "holes" in a regular DW 
structure {k = L/2- 2)0 

For a finite cyclic chain we determine the boundaries from equations: 

• FM-FP: e = 21 cos ^ 21 

• < 2 >-FP: e = -21 (cos ^ + cos ^) ^ -AT 

For an infinite chain we may introduce q = lim^^^oo k/L, the concentration of DW's, which 
may be varied from to 1/2. The analogue of (0) in the limit L ^ oo 

E r^F dk 

^ = qe-2{l-q)T — cosfc (14) 

L Jo IT 

where kp = vrg/(l — q). Differentiating over q in shows how q changes with e/pi: 

— = sin kp-\ cos kp . (15) 

P TT 1 — 

The first excited states whose energies are numerically calculated in Section III are non- 
trivial even in the framework of the free fermionic approach (cf Eqs.( |C. and ( |C.4| )). In 
Appendix C we derive Eq.( |C.2|) for the energy AE{Q) as a function of the wavevector Q. 
Such a fermion-hole excitation leaves the DW number unchanged. It becomes gapless at 
certain Q's when L oo. In this limit the only non-zero contribution to AE{Q) is given 
by Eq.O: 

{\sm(x - kp)\, 0<x<n/2 \0\ ( kiA 

' ^ ^' / ; a; = M 1 + ^ (16) 

|sin(x + MI, 7r/2<a;<7r ^ V ^> 

which is shown for kp <7r/2 in fig. 0. If kp >7r/2 one has to substitute it — kp for kp. The 
fermion-hole excitation plot is symmetric when a momentum transfer, Q, changes from to 
2kp{l + kp/n)^^. It formally differs of the results predicted by Eqs. (|C.3D and ( P.4|) . These 



occur due to a contribution of order 1/L in the first two terms in tlie r.h.s. of Eq.( |C.2| ). 
However, this contribution is still important when the finite-size systems are analyzed. 

This is noteworthy, that the contraction of the wavevector, 2kp 2kp{l + kp/Tr)^^, 
when a fermion is transfered from one Fermi-point to another, can be attributed to strong 
"anti" -correlated properties of our spinless "free" fermions, that forces for using a fictitious 
lattice description. 

As has been shown ini the FP should exhibit a power-like decay of correlation functions, 
e.g., 

{a- a (X cos Tcqx (17) 

with p=(l — g)^/2. This g-dependent p is also a result of the fictitious lattice contraction 
which increases with q. However, the FP cannot cover the whole range < q < 1/2: 
it appears to be unstable at smaller g's and transforms into the PM state through the 
mechanism of the Kosterlitz-Thouless transition at p ~ l/4i, corresponding to a wavevector 



gpM_FP = 1 - 1/V2 ^ 0.292 (18) 

Hamiltonian (^ allows to go beyond the first order in P which has been discussed above. 
A non-trivial contribution in the second order in P to perturbation theory results in ap- 
pearance the terms, which create (annihilate) DW's on next nearest sites. In fact, a couple 
of DW's on the sites, say, j and j + 2, can be created as a two-step process: First, if the 
high -energy excitation is virtually created by Tj^T^+ii then a many- fold degenerated ground 
state can be restored by 7y^_i'r,^2- The second possibility is in a sequential process: Tj^iTj^2 
and Tj^iTj^. Perturbation theory in the P^-order also extends the DW hopping terms up to 
next nearest neighbors, but this extension is out of relevance. 

Thus, in addition to Hamiltonian (P) we obtain 



^e? = -7 r;r;+2 + rtrU + r^rt + HC | (19) 



which can be also written in terms of fermions: 



^?ff = -7 c]cU + 4cLi + 44 + HC j (20) 

We put 7 = 2r^/l^ in Eqs.(|T9D-(pOD and HC denotes Hermitian conjugate. 

In spite of a smallness of Ti^g as compared to 7i|:g the former significantly influences on 
the critical properties of the FP: Tigg forces the FM-FP transition to be of the Ising type 
and trasforms this, in fact, into the FM-PM. This results in a gap opening in the PM phase. 
The excitation spectrum in a low DW density limit behaves as {k here is the wavevector) 



;e-2rcosA;)2 + (27sinA;)2 (21) 

which remains of a single-minimum kind within a narrow interval 2r < e < 2r - 27 VF, 
then it develops in a double-minimum curve. The gap value, in general, is oc F^ 

The situation on the FP-< 2 > boundary is different: This state is characterized by a reg- 
ular DW structure and elementary excitations driven by Ti),^ are the four - "leg" dislocations, 
which must be irrelevantlii. This "leg"- number can be easily illustrated in terms of the A 
and B particles. The FM vacuum is unstable at the FM-PM transition with respect to Tigg 
which transforms that vacuum state . . . BBBB . . . into the set of . . . B . . . BAAB . . . B . . . 
functions, that means p = 2. On the contrary, the < 2 > phase is described as the 

vacuum. According to the definition when A (or DW) disappears, it cre- 

(2) 

ates a couple BB. Hence, the elementary excitation due to Tigg can be represented as 
. . .A. . . ABBBBA . . .A. . ., which is irrelevant (p = 4). 

A fictitious lattice as it has been introduced ini is better realized in the case of infinite 
chain with no periodic boundary conditions imposed. The fictitious lattice sites are enumer- 
ated as m = m — ri{m), where ri{m) = J2i<m is the total number of fermions from the left 
of the m-th site of a real lattice. For fermions on the fictitious lattice we accept the notation 
cr (c~). TYgg is formally unchanged: 

= Y.^^^- r(4cf^+i + 4+1^-) (22) 

m 

while Ti^Q takes a form: 



= E + 1) + + 1)) (23) 

m 

where the operator P^{m) {P_^{rn)) translates all the fermions from the right of m by two 
sites to the right (left). 

III. EXACT DIAGONALIZATION OF FINITE SYSTEMS 

The starting point for our numerical investigation of the ground-state properties of the 
ANNNI-model (1) is the effective Hamiltonian Ti^g without creation/annihilation of domain 
walls (H), which we call the free model and the effective Hamiltonian + > which 



includes creation/annihilation of domain walls (19) and which we call the complete model. 
The Hilbert space for these models consists of all configurations of the original ANNNI- 
model, which obey the constraint discussed in the last section. This constraint reduces 
the dimension of the Hilbert space considerably so that much larger system sizes can be 
diagonalized (even if one uses the translational and spin flip symmetry of the original ANNNI 
Hamiltonian to block diagonalize it first). 



A. Methodology 

We reformulate the constraint in such a way that it becomes suitable for a numerical 
implementation. The domain walls can be identified with particles being able to hop to 
the left or right provided the constraint will not be violated by this move (moreover at 
most one particle can occupy a single site in the dual lattice). Since a domain wall at 
bond j comes always with bond j + 1 free of domain walls we call this combined object 
an ^-particle situated at bond j. If no domain wall is at bond j we call it a S-particle 
provided no other domain wall occurs at bond j — 1- In an obvious notation the following 
particle configuration, domain wall (i.e. r-) configuration and spin configuration (in the 
(Tz-representation) correspond to each other 
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. ITiiTTTTii ■■•) 

\AABBA---) = llOlOOOlO---) = ^"^^""^^ (24) 

liTTiiiiTT---) 

Moving a domain wall from bond i to i ± 1 means moving an ^-particle from site i to 
which is only possible, if at site i ± 1 is a i3-particle. 

T^T,-\BA---) = \AB---) 

(25) 

T{T+\AB---) = \BA---) 
Thus in the particle representaion the above mentioned constraint is already contained. 

(2) 

Analogous remarks hold for the creation/annihilation of domain walls described by Ti^g : in 
the particle formulation that means that ^-particles can be created in pairs in place of four 
consecutive i3-particles. For instance 

T+T+\BBBB ■ ■ ■) = \AA---) 

(26) 

T^TilAA---) = \BBBB---) 

Because of the periodic boundary conditions we have to discriminate between the cases with 
and without a domain wall at bond L. We denote the first group of states with a prime, for 
instance (in the case of L=8): 

\ABBAA)' = 101000101) (27) 

(note that in this notation the rightmost particle is always an ^-particle), whereas those 
without prime denote states without domain wall between L and 1, e.g.: 

\ABBAA) = 1 10001010) (28) 



This state is simply a circular left shift of the primed state in ([27|), in fact for each primed 
state there is a unique unprimed state that can be obtained from the former via a circular 
left shift. However, there are of course more primed states than unprimed ones. 

Thus we consider different primed and unprimed subspaces characterized by the number 
of ^-particles (i.e. number of domain walls), which is conserved under the action of Ti^g 
(note, however, that the latter mixes states of the primed and unprimed subspaces by moving 
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domain walls to or from the bond L. The dimension of these subspaces is simply given by 
the number of different possibilities to distribute ua and ua — 1 particles on L — and 
L — ua — i sites, respectively. 



dim„^ = 



dim' 



L-ua 

Ua 

L — ua — i 

UA-l 

The dimension of the whole Hilbert space we are considering is 

^ /L-nA\ /L-UA-l 
dimn^ff = 1 + 2^ < 

nA=l 



(29) 



(30) 



Ua ' ' ua-I 

which is a much smaller number than the dimension of the original Hilbert space, which is 
2^. 



L 




2^ 


dim^^ / 2^ 


8 


23 


256 


8.984375e-02 


12 


162 


4096 


3.955078e-02 


16 


1103 


65536 


1.683044e-02 


20 


7563 


1048576 


7.212639e-03 


24 


51842 


16777216 


3.090024e-03 


28 


355323 


268435456 


1.323681e-03 


32 


2435423 


4294967296 


5.67041 le-04 



Table I: Dimensionality of the reduced Hilbert space of our effective Hamiltonians (second 
column) and the dimensionality of the original Hilbert space (third column). 



From table I it becomes obvious that the storage requirements for diagonalizing the 
effective Hamiltonian is significantly smaller. This is still the case if one uses all symmetries 
of the original ANNNI Hamiltonian to blochdiagonalize it first (by which can reduce the 
storage requirement by roughly a factor 1/L). Note that in the FP-phase it is not a priori 
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clear in which wave number sector of the original Hilbert state the ground state is situated. 
Hence all of them have to be considered, as has been done inEl for lattice site up to L = 16. 
We were able to diagonalize easily L=32 systems on workstations with a reasonable RAM 
without reading or writing to the hard disk. 

In order to enumerate the states within the subspaces efficiently we recur to a scheme 
that is frequently used in the context of quantum spin chains confined to subspaces with 
constant magnetization^. Let us fix L to be a multiple of 4 (so that the periodic boundary 
conditions are fully compatible with the ground state in the anti-phase <2>). Let '?/'„^(n) 
be an unprimed state with ua ^-particles: 

V^„^(n) = \A{pi)---A{pnA)) = \B_^AB_^A ) , (31) 

where pi G {1, . . . , L — wa} denotes the position of the i-th ^-particle (counted from the 
left). Then the following definition yields a one-to-one correspondence between the possible 
configurations and a number n G {0, . . . , dim„^ ~ 1} • 

^ / Pi — 1 \ /■ i — \ \ 

n = Y,[ ) with ( j = . (32) 

The same definition is used for primed and unprimed states, which we denote by '?/'„^(n) 
and Tp'^^iji), respectively. The Lanczos routine we use to calculate the ground state and first 
excited state generates the Hamiltonian each time it is needed, for this reason we need to 
know how the various hopping, creation and annihilation operators act upon the basis states 
we have chosen. This can either be done by a hashing technique, which is used frequently for 
arbitrary quantum spin chains, or by explicitely calculating the number of the transformed 
state if possible. Fortunately in our case the latter is straightforward, and in Appendix D 
we list all relevant formulas that we need to generate the non-zero matrix elements for 
and T^eff • 

With the help of nowadays standard Lanczos routines we calculated the ground state and 
the first excited state of systems of size up to L=32. The effective Hamiltonians we derived 
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are expected to be good approximations to the original ANNNI model for small values of 
r/J. We confined ourselves to the values T/J = 0.05, 0.1, 0.2 and 0.5, where we calculated 
all physical quantities of interest for K-values in intervals of 5 ■ 10~^. We performed extensive 
checks of our code by comparing our numerical estimates for the ground state energy and 
the gap with the exactly known values for the free model (^. We also compared our results 
obtained for TY^g + 7^^^ at small values of F/J with those for the original ANNNI model 
and found no significant deviations. 



B. Results 

In order to check the quality of the free fermion description we first calculated the ground 
states of the original ANNNI model (|I])-(0) for a modest system size (up to L=20). The 
"classical" energy of a state 

£^class = imM) (33) 

which is -Eciass = —J'iJ2i SiSi+i — nY^i SiSi^2) for an eigenstate of the a^ operators contains 
the information on the average number of domain walls (or ^-particles) in the state ip, if 
this last does allow neither ... ||| . . . , nor . . . |t| . . . to appear. With such a constraint 
imposed on a state with exactly ua G {0, 2, 4, ... , L/2} domain walls, each of them costing 
an energy e = —(4k — 2) J' with respect to the FM state, we have 

E,i^JJ+{l-ti)L = nA-{4:ti-2) (34) 

Thus the comparison of the l.h.s. which we call DE, with the set of straight lines provides 
a measure of the average number of ^-particles plus an indication of the appropriateness of 
the free fermion concept in this context. In fig. we show the result for various values of F 
and we see that the smaller F the better the agreement of various parts of the Dii^-curves. 
Furthermore, for increasing F the ground state is more a superposition of various particle 
eigenstates in the vicinity of the FM-PM transition. 
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Now we turn our attention to the effective Hamiltonians. The average number of domain 
walls, simply given by (ra^) = {ip\ Tt'T^li') is calculated for ip's of the ground state and of 
the first excited state of the complete model (|^) + (|T9|) . For F not too large one observes well 
defined regions with constant value for ua in the ground state and the first excited state. 
See fig. for an example. The number of ^-particles increases monotonically (in steps of 
2) for increasing k. The points where the particle number jumps are of special interest. 
The particle number of the first excited state {ua}! jumps first abruptly (approaching the 
transition points from either side) and then it changes roles with the ground state. As a 
consequence the gap (i.e. the energy difference between ground state and first excited state) 
gets very small here. 

Note that in the FM-phase the gap does not vanish exponentially with system size for 
the effective Hamiltonians, because in the particle representation the two degenerated states 
with all spins up or all spins down are represented as one state. Therefore the gap closes 
in the infinite system only at the FM-PM transition. On the other hand the gap stays 
zero (or exponentially small for finite sizes) throughout the < 2 > phase also for the effective 
Hamiltonians: The 4-fold degeneracy there is only reduced by a factor two via the elimination 
of the spin-flip symmetry and a degeneracy between corresponding primed and unprimed 
states is left. 

As a significant difference of the complete model with respect to the free model we 
note that at the special K-values, where the particle number changes and which can be 
calculated exactly via the formula ([T3|), the gap of the free model (§) closes completely, i.e. 
^E^nT_^nA±2 — 0- complete model the gap-value on the boundary, say, — + 2, 

can be easily estimated if we confine our consideration with these two competing states only, 
that results in 

= 2(^„,+2|HSVnJ (35) 

These special gap- values increase with F because of the F-dependent T^es • particular 
more pronounced is the gap increasing for the lower K-values. In this range the boundary 
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Ha ^ ua + 2 cannot be considered as well isolated from other "neighboring" states, i.e., 
UA — 'i, Ha + 4:, etc. This observation will turn out as a hint for the existence of the PM-FP 
transition. For higher values of P/J the resulting picture is therefore slightly different: the 
transitions smear out and {nA){K) and AE become smoother as can be seen in fig. ^. 

Let us interpret this picture cum grano salis: One might tentatively locate the FM-PM 
transition (for P/J = 1/5) roughly at k ~ 0.41, where the gap should approach zero like 
AE ~ 1/L (since the transition is expected to be in the Ising universality class, where 
u = z = 1). Between k = 0.41 and let us say k, ^ 0.48 the average particle number increases 
from zero to 8 (for L=24), but the individual transitions observable in fig. ^ melt together 
to form a rugged plateau. Only when ua gets larger than some value (whose significance we 
will clarify later) the gaps at the individual transitions try to close again. With increasing 
system size these gaps (for k, larger than roughly 0.48) melt together, too, but in the limit 
L — s> oo they will form a curve A(k) = in this region, which is simply the gapless floating 
phase with a quasi-long range, i.e., algebraically decaying spin correlations. 

It is obvious that in order to observe this scenario in its pure form one has to go to 
enormous system sizes, which is not feasible yet. Nevertheless we can clearly demonstrate 
the qualitative difference between the PM phase and the FP phase by exphcitly studying 
the spin correlations in both regions of the phase diagram for intermediate system sizes. 

The spin correlation function is defined via 

Cir) = {i^\atat^,\^) 

L/2 dim„^ ^ L (36) 

= I] ^ V'^^(n) ■ -^^i[^„^(n)] ^i+^V'nAH] + primed states 

n^=0 71=0 i=l 

Here i5['?/'n^(n)] (5 = Si, ... , Sl) means the spin configuration that is equivalent to the state 
number n with ua particles (since there are always two of them we choose the one with the 
first spin up 5*1 = +1). Of course we have to initialize such a mapping in our program once, 
afterwards this table can be used whenever correlations have to be calculated. 

First we take a look at the structure function since this directly relates to the particle 
number ua discussed above. We define it as follows: via 
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/.(«^) = COS (g(^ - j)) {M^t^'M) (37) 
for q = n/L (n = 0, 1, . . . , L — 1) so that 

C(r) = ^cos(gr)/, (38) 
g 

In fig. D we show the result for /^(k) for T/J = 1/20 in comparative plots for the free and 
the complete model. For fixed wave number q the structure function /^(k) is simply a step 
function for the free model, the plateaus located at the K-intervals with constant particle 
number n^. For fixed k the wave number g^ax with the maximum amplitude fq^^{i^) > fq{i^) 
for q 7^ ginax is related to the particle number via 

gmax = ua/^L (39) 

(note that Ua is a good quantum number in the free model). For small F/ J we observe that 
the steps get rounded, but nothing dramatic (for these system sizes) happens. If we increase 
F/J, as is done in fig. ^, the neighboring plateaux begin to mix (to melt, see above). They 
also shift their location to larger values of k. However, these changes seem to become less 
significant as soon as gmax is larger than 1/4 (actually, as stated in the last section in Eq.(^), 
the theoretical estimate for this relevance-threshold is approximately q = 0.292 > 1/4). 

This statement finds its strongest support by looking at the spin correlation function 
(^), shown in fig. ^ and ^ for L=32, directly. For the free model of finite size L with 
periodic boundary conditions one would according to (|1^) expect that 

G{r) = a cos{rnq^,^) ■ {r-P + {L - r)-P} with p = ^(1 - gmax)' , (40) 

where a is a fit parameter and gmax has to be determined from the structure function. In 
order to resolve the correlations over as large as possible distances we took here the largest 
possible system sizes. For L=32, however, we had to confine ourselves to a smaller number 
of parameter values. In the case F/J = 1/20, in which the differences between the complete 
and the free model are not too large, we took simply the middle of the plateaus of the 
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structure function, whereas for F/J we took the location of the maxima of fq{n) shown in 
fig- I- 

What is shown in fig. ^ and ^ is a comparative plot of C(r) and G{r) with a = 



J2r C'^i.'f)/ I]r G'^i.f)^ Cl{r) = G{r)/a. We observe that for gmax > 1/4 the free model corre- 
lation function G{r) fits C(r) in an excellent way. For gmax < 1/4, however, one recognizes 
sigificant differences between G{r) and C(r), most dramatic for the smallest wave numbers, 
which mean closest to the FM-PM transition of the complete model. The gmax = 1/12 and 
Qmax = 1/6 curves for F/J = 1/5 definitely decay faster than algebraic. We obtained an ex- 



cellent fit (shown in fig. |ToD for the Qmsix = 1/16 curve by a superpostion of an exponentially 



damped g = and q = 1/16 oscillation: 

a=o.4085;r/j=5(r) = (1 - a + acos(r7r/16)) ■ (e"^/? + g-^^'^-r-)/?) (4^) 

with a = 0.21 and ^ = 3.6. For us this is ample evidence that for fixed F/J one enters 
first a paramagnetic phase with exponentially decaying correlations by increasing k from 
the FM-phase. Only when q^ax gets larger then the above mentioned value, one enters the 
FP-phase, which for the free model extends over the whole region between the < 2 > and 
FM phases. 



IV. SUMMARY 

Below we summarize a few important points of this work. 

• Instead of the original Id ANNNI model in a transverse field we consider a reduced 
model which we show to be a very reasonable modification when the competition 
parameter |k — 1/2| as well as the quantum parameter T/J' are small. 

• The effective Hamiltonian Ti^g + Ti^g is most easily visualized with the A — B particle 
representation. The latter has been adapted to the periodical boudary conditions and 
efficiently used in numerical calculations. 
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• Analytical calculations here were done not for re-discovering aftera the essential 
physics, accompanying the scheme of phase transitions FM-PM-FP-<2> in infinitely 
long chains. Most analytical results are derived for further applications in a numerical 
scheme. 

• The exact numerical diagonalization technique (Lanzcos algorithm) was used on the 
restricted basis states, all of them have been enumerated, as well as all the non-zero 
matrix elements generated by Hamiltonian of the complete model were stored. Along 
this line we perform calculations for systems of size up to L=32 with standard work- 
stations. 

• The set of Figs. 5-10 convincingly illustrates a different origin of the critical behavior 
at larger and smaller k's. Most likely, the PM-FP transition takes place close to or 
even at a modulation q given by the theory of the standard Kosterlitz-Thouless-like 
transition. 
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Appendix A 

We treat domain wall operator, defined on the lattice site j of the dual lat- 

tice. It takes two values, ±1: +1 (—1) signals of absence (presence) of a domain wall. 
Mathematically, this can be expressed by 

Because of the periodic boundary conditions we also define 

n' = ^V2^^3/2 and tI = al/^al_y^. (A.2) 

Also evidently, that 

= ^1-1/2^^1+3/2- (A-3) 

For x-components we accept the following definitions: 

^J+i/2 = (A-4) 

and 

'^i/2 = rirl and (Tl_y2 = t^_^t2. (A.5) 



With using Eqs.( |A.l| ),( |A.3| ) and ( |A.4| ) we can easily obtain the form of Eq.(^ from (|I]). 



This transformation set allows to calculate the correlation functions in terms of r's, for 
example: 

«-i/2<+i/2) = ( n (A-6) 

j=ri 

Appendix B 

In Eg. (pup m's are supposed to be arranged in order, satisfying the following constraints: 
nil < ^2 — 1, ""^2 < "^3 — 1, . . . ,mk-i < rrik — 1, rrifc < m\ + L — 1. We transfer these 
constraints on the amplitudes /: 
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/(. ..,m,m + l...) =0 (B.l) 

/(1,...,L) = 0. (B.2) 

In spite of the irregular 1 ^ L hopping term in Eq.(^ the eigen amphtudes / satisfy the 
regular equations: 

(E - ke)f{mi,m2, rrik) = -'^ (/("^i + a, ms, . . . , m^) + 

a=±l 

f {mi, 1712 + a, . . . , mfc) + . . . + /(mi, ms, . . . , + a)). (B.3) 
Using a Bethe substitution for /'s (see Eq.([TT|)) we obtain 

k 

E = ke-2Tj2(^osqj. (B.4) 



Eqs.QBlD and (^) yield 

LAj^.e'"^ +U^..e''' =0 (B.5) 

and 

respectively. In turn, from Eqs.( |B.5D and ( p.6|) one can arrive to 

27r , , 

Qi - Qj = j;:ri^i3 y^-^) 

where riij are integer numbers. Additional equations imposed on g's may be obtained from 
Eq.( [B.3| ) at m = 1. Formally, it is equivalent to /(O, m2, . . . , m^) = /(m2, . . . , m^., L), that 
is 

^^,s...t = Os...t.e*'"'' (B.8) 
which after a simple algebra results in 

odd 1 k 

l^oddj^ the set of odd integer numbers. 

The content of this Appendix may be summarized in this last Eq. (|B.9|) and the form of 
the the eigenfunctions ([T^). 



Appendix C 

For the ground state of k DW's the set {rf'^'^} in Eq.( [B.9| ) counts all the equidistant odd 
numbers from —k + liok — 1. Let us take away one number from that set, say, n"^^ = 2z/_— 1, 
with —k/2 + 1 <u_ <k/2 — 1, and add another number 

n°dd= 2u+- 1, with u+>k/2. With 
this definition ^i'^'^ = 2(z^+ — z^-) and after summation in both parts of ( [B.9| ) we get 

271 (u, — u) 



(C.l) 



The wavevectors in Eq.( [B.9| ) are now determined and the energy of such an excitation 
counted from the ground state energy takes a form: 

A^(Q) = -2 5:cosg,-Eg3 = -2 E cos ^^'^ ~ ~ ^ - 

i v=-kl2+l ^ 

{2u^-l)n-Q , ^ {2u^-l)n-Q , ^ ^ 7r(2z/ - 1) 

-2 cos ; -, h 2 COS ; ; h 2 ^ COS ■ 



L-k 



L-k 



v=-k/2+l 



, Q \ nk 

2 1 — cos - — - sin 



—2 sin 



Q 



-2 cos 



L-k 
Q 



sm 



L-k 
7r(2z/+- 



L-k 



L-k 

I TT 

sm 



L-k 



L-k 



sm 



7r(2z/_- 1) 



L-k 



7r(2z/+-l) 

cos ; COS 

L-k 



L-k 
7r(2z/_- 1 
L-k 



(C.2) 



Using Eq. (|C.2|) we can analyse the lower lying excitations 



1. n°'i'^= k-\, nX'^= + l ^Q = 2tt/L. 

The excitation energy with a vanishing momentum transfer Q in the L oo limit 
reads in the leading order in 1/L (this contribution comes from the last term of the 
r.h.s. of Eg. (1021) ): 



A 4:71 . Tik ^ A . tt/c 
A£'((5) ^ sm - — - = 2Agsm 



L-k L-k 



L-k 



(C.3) 



where Aq is a g spacing (cf Eq. (|B.7|) ) 



2. ^odd^ ^ Q ^27Tk/L. 

All the terms in the r.h.s. of Eq.( |B.7D contribute in the leading order in 1/L. This 
excitation occurs at the energy 
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AE{Q) « 2A, (l - 1^ 



(0.4) 



which differs from ( |U.ci| ) . 



3. For not a special value of Q a leading contribution, 0(1), comes from the last term of 
the r.h.s. of Eg. ( |C.2| )). To leading order it can be written as: 

AE{Q) =Asmxsm{x + 2TTiy^/{L-k)), x = (C.5) 



We used (|C.1|) to obtain the form of Eq. (|C.5|) . Then in order to select the lowest 



energy values at fixed Q we consider following inequalities: 

• X < kp = 7[k/{L — k). The possible range of z/_'s becomes LQ/{2'k) — k/2 < 

<k/2. The minimum of (|C.5|) is reached at the lower limit, that results in 
AE{Q) = 4 sin x sm{kF — x) 

• kp <x <7r — kp- Now a possible range of z/_'s is —k/2 < z/_ < k/2. Two extreme 
possibilities should be checked, one arrives to the form AE{Q) = 4sina;sin(x — 
kp), another results in AE{Q) = 4sina;sin(A;ir + x). The former realizes the 
minimum at kp <x <7t/2, the latter is correct at 7r/2<x<7r — kp. 

• TC — kpKxKTT — > —k/2<u^<L — 3k/2 — LQ/(2'n-). A true minimum corresponds 
to AE{Q) = 4 sinx sin(— a; — kp) 
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Appendix D 

In this Appendix we list the formulas that determine the action of various hopping, 
creation and annihilation operators on states ■?/'„^(n) and ip'^^{n'). First remember the 
definition of our notation (0). In what follows the numbers n and n' are always given by 

/ p,- - 1 \ "A-i / p,- - 1 ^ 

n = E and n'= Y.[ ) , (D-1) 

First we consider the hopping term occuring in (|]): 
Hopping to the right t.+ 

The matrix elements of the operator X^i^i t^t^+i in the particle representation is non-zero 

whenever it is possible to move an ^-particle to the right. This means there has to be a 
i G {1, . . . in such a way that pj^i > pj + 1, compare with (|25D. We have to take 
special care of the case i = L — 1 or L (i.e. hopping to or from the periodic boundary), when 
an umprimed state transforms into a primed state and vice versa. 

i e {1, . . . ,L - 2}, i.e. < L - ha- 

tlJnA'^) = I ■ ■■A{p,)B-- ■) ^„Jm) = I ■ --BAip, + 1) ■ ■ ■) 

fPj-l\ (D-2) 



m = n + 



J -I 



The same for a primed state ipn^{n'). 
i = L-1, i.e. pnA = L- ua: 

^ljnAn) = \B---A) — . ^;jm') = |---S^)' 



m 



3=1 J 



i = L, i.e. all primed states: 



^A^^fPj-l\ (D-4) 
m = 2^ ( 

j=i ^ J + 1 
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Hopping to the left Tj^r^^i 

i e {1, . . . ,L - 2}, i.e. < L - ua: 

^IJnAn) = I ■ --BAip,) ■■■) V^„^(m) = I ■ --Aip, ~1)B-- •) 

J-1 



m = n . 



The same for a primed state ipn^i'i^-') ■ 
i = L, i.e. pi = 1: 



i=2 ^ J - 1 



z = L — 1, i.e. aU primed staes: 

ij'^^in') = \ ■ ■ ■ BAY — ^nAm) = \B---A) 



m 



L~nA-l\ ''A^^^Pi\ (D.7) 



Next we consider the creation and annihilation operators occuring in the complete model 
via ([T9|). Again we have to take special care of the cases in which a domain wall at the bond 
linking site 1 and site L is created or annihilated. 



Creation of domain walls t{~t^2 



The non-zero matrix elements of the operator X^i^i '^i''^i+2 have to be determined via the 



rule (|26|) . Let us denote with py a position between two succesive ^-particles at positions 
Px and px+i (with p^+i — > 4) where a new pair of .4-particles can be created. 

Py e {1,. . . ,nA- 3}: 



25 



V„.,(n.) = \A{p,) ■ ■ -^(p,) \B{p,)B{p, + l)B(p„ + 2)B{p, + 3)] ^(p,.+,) • --Aip^J) 
V-n^ (m) = |^(pi) • ■ ■ ^(Pi) ^(p„)^(p, + 1) A(j},+i-2)---A(j}n^-2)) 



j=i ^ j ^ ^ X + 1 ^ ^ X + 2 ^ j=x+i ^ j + 2 



I ^ 



The same for a primed state ip'^^ (n') with replaced by — 1 in the last sum. 
^ = - 1) • --Aipn, - 1)A{L -HA- ?>)A)' 

j=l^ j ^ ^ UA + l 



Py = L - UA-l: 



Py^ L- Ha- 



^^,{n) = \BBA{p^)---A{pnJBB) 

i'nAim) = \A A{pi - 1) ■ ■ ■ A{pnA - 1) A) 

^ = ) + ( 

3=1^ j + 1 ^ ^ nA + 2 



^„^(n) = \BBBAip,)---A{pnA)B) 

= 1-4(1) Aipi - 2) . .-Aipn, - 2) AY 



m' = 

j=i ^ J + 1 



(Di 



(D.9) 



(D.IO) 



(D.ll) 



Annihilation of domain walls t^_^_2 

Now let there be two successive ^-particles at position p^ and Px+i = Px + I, so that they 
can be annihilated by t^t^^2 fo^ some suitably chosen i. 
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Px e {1, . . . ,L - ua - ^Y- 

^„^(n) = \Aipi) ■ ■ ■ Aipx-i) Aipx)Aipx + I) Aipx+2) ■ ■ ■ Aipn^)) 

(m) = \A{pi) ■ ■ ■ AiPx^i) B{px)B{px + l)B{px + 2)B(p, + 3) A{px+2 + 2) • • • + 2)) 

^ / - 1 \ "4 / Pi + 1 



j=l ^ j ^ j=x+2 ^ j -2 



The same for a primed state (n') with replaced by — 1 in the last sum. 
Pj; = L — — 2 for a primed state: 

<Jn') = \A{p^) ■ ■ ■ A{pnA-2) AAY 

il^nAm) = \BA{pi + l)---A{pnA-2 + l)BBB) 

nA-2 



(D.12) 



(D.13) 



.7=1 ^ T ^ 



m = 



Px = L — Ha — i ior an unprimed state: 

^ V'^^lm) = \BB A{p2 + 1) • --AiPnA-i + 1) -BiS) 



n^l Pj 
i=2 



(D.14) 



Px = L — ua for a primed state: 

^ ^,^(m) = |i3S-B^(p2 + 2)---^(p,^_i + 2)-B) 



(D.15) 



i=2 ^ j - 1 



In order to determine the non-zero matrix elements for creation/annihilation it is thus nec- 
essary to scan the whole particle configuration corresponding to 'il){n) and to record the 
allowed values for y and Px and Py. An alternative technique would be to try to store the 

(2\ 

(very sparse) matrix for . However, as in most cases, not the computational speed but 
the storage requirement is the limiting factor for the maximum possible system size. 
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FIGURES 



FIG. 1. The one particle excitation energy AE{Q) for the free model from (16) versus the 
contracted wave vector Q ■ (1 + kp/7r)/2 for different values of kp = — q), with q the domain 

wall density (which can be evaluated as a function of k and F/J via equation ([Tsl)). kp = 
corresponds to the FM-phase, kp = 7r/2 corresponds to q = 1/3, and for kp between tt/2 and tt 
(q = 1/2, i.e. <2>-phase) the evolution is simply reversed. 

FIG. 2. Results of exact diagonalization of the original ANNNI model (|^)-(|2|) for L=16 and 
various values of F/J: Expectation value of the qunatity Edass/ J + (1 ~ ^^)L defined in the text. 

FIG. 3. Results of exact diagonalization of the complete model (||t)+(|^)- Shown are the 
results for the expectation value of the particle number in the ground state and the first excited 
state together with the gap A^;. It is L=24 and F/J = 1/20. 

FIG. 4. The same as in fig. | but with F/J = 1/5. 



FIG. 5. Comparision of the structure function fq{n), eq. (|37|), for the free and the complete 
model. It is L=24 and T/J = 1/20. Note the different scale on the y-axis in the plot for q = 7/12. 

FIG. 6. The same as in fig. | but with F/J = 1/5. 

FIG. 7. Comparision of the correlation function C(r), eq. (^), of the complete model with 
G(r) given in eq. (|^). It is L=32 and F/J = 1/20. 

FIG. 8. The same as in fig. but with F/J = 1/5. 

FIG. 9. Structure function /^(k), eq. (^), for the complete model in the case L=32 and 
F/J = 1/5. 

FIG. 10. Fit of (|l|) to the correlation function C(r) for k = 0.4085, the maximum of the 
fq=i/i2 curve shown in fig. It is L=32 and F/J = 1/5. The fit parameters are a = 0.21, ^ = 3.6. 
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ANNNI chain: J"^ = 0.1; T = 1; L = 16 
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